A Measurement of the Coulomb Dissociation of 8 B 
at 254 MeV/nucleon and the 8 B Solar Neutrino Flux 
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We have measured the Coulomb dissociation of B into 
7 Be and proton at 254 MeV/nucleon using a large-acceptance 
focusing spectrometer. The astrophysical S17 factor for the 
7 Be(p,7) 8 B reaction at E c . m . = 0.25 - 2.78 MeV is deduced 
yielding Sir(0) = 20.6±1.2 (exp.) ± 1.0 (theo.) eV-b. This re- 
sult agrees with the presently adopted zero-energy Si 7 factor 
obtained in direct-reaction measurements and with the results 
of other Coulomb-dissociation studies performed at 46.5 and 
51.2 MeV/nucleon. 

PACs: 25.40.Lw, 25.60.-t, 25.70.De, 26.65.+t 



The precise knowledge of the solar thermonuclear fu- 
sion of 8 B (from 7 Be plus proton) is crucial for estimat- 
ing the 8 B solar neutrino flux and the predicted neutrino 
rates in terrestrial neutrino measurements. The relevant 
7 Be(p,7) 8 B cross section a(E) is parameterized in terms 
of the astrophysical factor Sn(E) which is defined by 
S 17 (E) = a(E)Eexp[2Trr](E)] where rj{E) = Z Y Z 2 % 2 /hv 
is the Sommerfeld parameter. The flux of 8 B solar neu- 
trinos is particularly important for the results of the 
Homestake, Super Kamiokande, and SNO experiments 
JIJ which measure high-energy solar neutrinos mainly or 
solely from the 8 B decay. 

Unfortunately, this cross section has not been known 
with sufficient accuracy for a long time, despite the fact 
that several comprehensive direct measurements were re- 
ported for the 7 Be(p,7) 8 B reaction k|. The main dif- 
ficulty in such experiments is the determination of the 
effective target thickness of the radioactive 7 Be target. 
This difficulty is reflected in the fact that the results 
of these measurements can be grouped into two distinct 
data sets which agree in their energy dependence but dis- 
agree in their absolute normalization by about 30%. In 
view of this discrepancy, experimental studies with dif- 
ferent methods are highly desirable. 



As an alternative approach one can measure the in- 
verse process, the Coulomb dissociation (CD) of 8 B into 
7 Be and proton Q . The CD yields are enhanced because 
thicker targets can be used and a larger phase space is 
available for CD. This method uses stable targets and 
thus is free from the difficulty of determining the effec- 
tive target thickness. On the other hand, direct (p/y) and 
Coulomb dissociation measurements have different sensi- 
tivities to the multipole composition of the photon fields. 
The E2 amplitude is enhanced in CD due to the large 
flux of E2 virtual photons, whereas it can be neglected 
in the (p,7) reaction. 

Recently, Motobayashi et al. have performed a CD ex- 
periment at E( 8 B) = 46.5 MeV/nucleon, yielding values 
for Sn in the energy range 0.6—1.7 MeV The ex- 
tracted (p,7) cross section is consistent with the results 
from the lower group of direct-reaction data points [Q-^) . 
Another measurement at 51.9 MeV/nucleon by the same 
group with improved accuracy led essentially to the same 
conclusion ||. 

In this article, we report on an experiment of the CD 
of 8 B at a higher energy of 254 MeV/nucleon performed 
at the SIS facility at GSI, Darmstadt, Germany. The 
present incident energy has several advantages compared 
to those used in Refs. ^||: (i) due to the strong for- 
ward focusing of the reaction products the magnetic spec- 
trometer KaoS |[o| was used for a kinematically complete 
measurement with high detection efficiency over a wider 
range of the p- 7 Be relative energy; (ii) effects that ob- 
scure the dominant contribution of El multipolarity to 
the CD, such as E2 admixtures and higher-order contri- 
butions, are reduced liljji^l ; (iii) the Ml resonance peak 
at E rc \ = 0.63 MeV is excited stronger and therefore can 
be used to check the accuracy of the invariant-mass cal- 
culation. 

A 8 B beam was produced by fragmentation of a 
350 MeV/nucleon 12 C beam from the SIS synchrotron 
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impinging on an 8.01 g/cm 2 beryllium target. The 
beam was isotopically separated in the fragment sepa- 
rator (FRS) jl3| and transported to the standard tar- 
get position of the spectrometer KaoS |l(J . The average 
beam energy of 8 B at the center of the breakup target was 
254.0 MeV/nucleon, a typical 8 B intensity was 10 4 / spill 
with 7s extraction time. Beam particles were identified 
event by event with the TOF-Ai? method by using a 
plastic scintillator with a thickness of 5 mm placed 68 m 
upstream from the target and a large-area scintillator 
wall placed close to the focal plane of KaoS. About 20 
% of the beam particles were 7 Be, which could however 
unambiguously be discriminated from breakup 7 Be par- 
ticles by their time of flight. 

An enriched 208 Pb target with a thickness of 199.7 (± 
0.2) mg/cm 2 and with an effective area of 20 x 22 mm 2 
was placed at the standard target position of KaoS. The 
reaction products, 7 Be and proton, were analyzed by 
the spectrometer which has a large momentum accep- 
tance of Pmax/Pmin ~ 2 and an angular acceptance of 
140 and 280 mrad in horizontal and vertical directions, 
respectively. Two pairs of silicon micro-strip detectors, 
installed at about 14 and 31 cm downstream from the 
target, measured the x- and y-positions and hence the 
scattering angles of the reaction products in front of the 
KaoS magnets. Each strip detector had a thickness of 
300 /im, an active area of 56 x 56 mm 2 , and a strip 
pitch of 0.1 mm. The energy deposited on each strip was 
recorded using an analog-multiplexing technique. This 
configuration enabled us to measure opening angles of 
the reaction products with a la-accuracy of 4.8 mrad 
which is mainly caused by angular straggling in the rel- 
atively thick Pb target. By reconstructing the vertex at 
the target with an ler-accuracy of 0.3 mm, background 
events produced in the target frame or in the strip de- 
tectors could be largely eliminated. The remaining back- 
ground events were found without the target under the 
same condition to be less than 0.5% of the true events. 

Momenta of the reaction products were analyzed by 
trajectory reconstruction using position information from 
the micro-strip detectors and two two-dimensional multi- 
wire proportional chambers (MWPC) which detected the 
protons or the Be ions close to the focal plane of KaoS 
with a position resolution of about 1 mm. 



A large-area (180 x 40 cm 2 ) scintillator wall |10| con- 
sisting of 30 plastic scintillator paddles with a thickness 
of 2 cm each was placed just behind the MWPC. It served 
as a trigger detector for the data acquisition system and 
as a stop detector for time-of-flight (TOF) measurements. 
To distinguish breakup events from non-interacting beam 
particles, the TOF wall was subdivided into two equal- 
size sections. The breakup events were characterized by 
coincident hits in the low-momentum (proton) and high- 
momentum ( 7 Be etc.) sections of the wall. The beam 
normalization was done by analyzing 8 B single hits on 
the high-momentum section (downscaled by a factor of 
1000) under the condition (from the tracking informa- 
tion) that they originated from the active target area. 



From the measured angles and momenta of the breakup 
products the p- 7 Be relative energy and the scattering an- 
gle 6s of the center-of-mass of proton and 7 Be (i.e. the 
excited 8 B) with respect to the incoming beam were re- 
constructed. 

The p- 7 Be coincidence yield is shown in Fig. [I] (a) as a 
function of the p- 7 Be relative energy. It is obvious from 
the figure that, contrary to the direct measurements, the 
Ml resonance at E rc \ = 0.63 MeV is not very pronounced. 
This is related to the energy resolution of our experiment 
(see below), and also to the relatively low sensitivity of 
CD to Ml transitions. 

To evaluate the response of the detector system, 
Monte-Carlo simulations were performed using the code 
GEANT. The simulation took into account e.g. the finite 
size of the Si strip and MWPC detectors as well as our 
inability to discriminate proton and Be for very small 
opening angle where both particles hit the same strip. 
Events were generated with probabilities proportional to 
the CD cross section calculated with a semi-classical for- 
mula |l||. For the Ml resonance at £ re i = 0.63 MeV, the 
(p,7) cross section calculated from the total and gamma 
widths measured by Filippone et al. @] was used. The 
non-resonant contribution was obtained by normalizing 
the El (p,7) cross section calculated by Bertulani pl| 
with a scale factor of 1.20. 

Further corrections in the simulation are due to the 
feeding of the excited state at 429 keV in 7 Be. We used 
the result by Kikuchi et al. ]l6| who measured the 7- 
decay in coincidence with the CD of 8 B. 

The histograms in Fig. [l](a) show the simulated 
El+Ml yields (solid) and El yields (dashed). As seen in 
this figure, the shape and magnitude of the experimental 
energy dependence are well reproduced. This indicates 
that the CD yield is well described by a combination of 
the Ml resonance and the pure El continuum. In the 
lower part of Fig. [l] we show that the total efficiency cal- 
culated by the GEANT simulation is high over the entire 
E le \ range covered in our study. From the Monte Carlo 
simulation we also estimate our relative-energy resolution 
to be e.g. a(E ml ) = 0.11 and 0.22 MeV at E mi = 0.6 and 
1.8 MeV, respectively. 

In order to estimate upper limits for a possible E2 con- 
tribution to our yields, we have analyzed the 9g distri- 
butions. Since we did not measure the incident angle 
of 8 B at the target, the experimental angular distribu- 
tions represent the 6$ distributions folded with the an- 
gular spread of the incident beam. In Fig. |^, we plot 
the experimental yield against the scattering angle 9g 
for three relative-energy bins indicated in the figure; the 
contribution from the measured beam spread is shown 
in Fig. H (d). The full histograms represent the results 
of a Monte-Carlo simulation for El excitation using the 
theory of Bertulani jl5| , normalized to the experimental 
yield. The Ml transition also contributes to the 0.5—0.7 
MeV bin with the same angular dependence. The sim- 
ulated E2 angular distributions are also shown by the 
dashed histogram. The angular resolution was estimated 
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to be 0.35 degrees (Iff), smaller than the observed widths. 
As seen in the figure the experimental distributions are 
well reproduced by the simulation using El and Ml mul- 
tipolarities, and no room is left to add an additional E2 
component. The corresponding 3ct upper limits for the 
E2/E1 transition amplitude ratio are calculated from our 
model as S B2 /Sei = 0.06 x 1CT 4 , 0.3 x 10~ 4 and 0.6 x 10~ 4 
for E ld = 0.3 - 0.5, 0.5-0.7, and 1.0-1.2 MeV, respec- 
tively. 

Recently, Esbensen and Bertsch have pointed out that 
interference between El and E2 components could play 
an important role in CD of 8 B and should manifest itself 
in an asymmetry of the longitudinal-momentum distri- 
bution of 7 Bc ]l7| . Analyzing measured 7 Be momentum 
spectra in term of this model, Davids et al. Jl8| have de- 
duced S E2 /S E i = 6.7±f;| x 10" 4 at 0.63 MeV which is 
much higher than our upper limits above. One should 
keep in mind, however, that such a strong E2 component 
need not necessarily affect the extraction of S17: the dy- 
namical calculation of Esbensen and Bertsch |l7| at 45.6 
A MeV shows that its effect can be much smaller than 
what is obtained by simply adding El and E2 compo- 
nents. Further experimental and theoretical studies are 
clearly desirable to clarify the magnitude of the E2 com- 
ponent in CD. 

From the data and the simulation presented in Fig. [l] 
(a) we can deduce the astrophysical Si 7 factors for El 
multipolarity by scaling the theoretical S17 factors by the 
ratio of observed and simulated counts. The resulting Si 7 
distribution is plotted in Fig. || as a function of the p- 7 Be 
relative energy, together with existing (p,7) and CD 
H results. The (p,7) results were scaled to a peak cross 
section of the 7 Li(d,p) 8 Li reaction of 147 ± 11 mb 
The binning of our data was chosen to be approximately 
equal to the FWHM of the E Te \ resolution. The errors 
shown in Fig. 3 result from summing in quadrature sta- 
tistical errors and those resulting from uncertainties in 
the momentum calibration, in the angular cutoff of pure 
Coulomb processes, and in the feeding of the excited state 
in 7 Be. Since we are mostly interested in the El com- 
ponent, the simulated Ml component of the E IC \ = 0.63 
MeV resonance was subtracted from the data. At low 
E Ic i, our results agree well with the lower group of the 
(p,7) results [|-| and the CD result of Kikuchi et al. [§. 
We observe some discrepancies at higher E le \, however, 
as shown in Fig. |3|. While we cannot explain the discrep- 
ancy with the RIKEN data [|l6| , since all of the assump- 
tions underlying their and our analysis are identical, the 
discrepancy to the direct measurements may be due to 
our neglect of an E2 component which, in a proper theo- 
retical treatment, might affect the high-i? ro i data points 
while leaving the low E re i data virtually unchanged (see 
right-hand panel of Fig. 11 in Ref. Jj"7|). 

Several theoretical models for 8 B have been proposed 
to predict the shape and magnitude of the S17 en- 
ergy dependence, such as one-body potential models 
(see e.g. |[l],[l5) for recent calculations) or more complex 
many-body (cluster) models (e.g. p0|-p3|). Our data in 



Fig. p| f ollow very closely the results obtained by Bertu- 
lani p5| over the entire range of energies as indicated by 
the solid curve which was normalized to our data points 
using a scaling factor of 1.20. The shape of the distri- 
bution predicted by the cluster model calculation p0[ , 
which was favoured by the recent (p,7) experiment of 
Hammache et al. Q, does not agree equally well with 
our results at the higher energies (dashed line in Fig. |^). 

To obtain the zero-energy astrophysical Si 7 factor, we 
follow Jennings et al. (2^] in fitting the theoretical en- 
ergy dependence only to data points below E le \ « 0.45 
MeV (see insert in Fig. ||) since this region should be 
largely free from uncertainties concerning nuclear excita- 
tions. From our two low E te \ data points we extracted 
Si 7 (0) = 20.6 ± 1.2 ± 1.0 eV-b where the first contribu- 
tion results from fitting the data points of Fig. |^ and the 
second one is related to the uncertainty in extrapolating 
to zero E re \. This result is compatible with the value of 
Si 7 (0) = 19.0±1.0±0.2 eV-b obtained by Jennings et al. 
p4j when fitting the combined data points of Refs. 
and also with the adopted value of Si7(0) = 19^ eV-b 
from Ref. |19[|. We note, however, that the question of 
possible backscattering losses of recoiling nuclei in direct 
(p,7) and (d,p) reactions |2[||26| is not yet solved. 

We conclude that we have demonstrated that high- 
energy Coulomb dissociation is very useful for determin- 
ing the astrophysical S-factor of the 7 Be(p,7) 8 B reaction 
at low energies. We have used an analytical formula from 
the literature [£4| to extrapolate our two low-energy S- 
factors to Si7(0) and obtain a value that is consistent 
with the most recent compilation p9[ . This supports 
essentially the standard-model prediction for the high- 
energy solar neutrino flux |FJ . 
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FIG. 2. (a, b,c) Yields of breakup events plotted against 8g, the 
scattering angle of the excited 8 B, for three relative energy bins. 
The full histograms show the results of a simulation taking into ac- 
count the measured angular spread of the incident 8 B beam (shown 
in d) and assuming El+Ml multipolarity. The dashed histograms 
show the simulated results for E2 contribution according to the 
calculations of Bertulani [ ^jf . 
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FIG. 1. (a) Yields of breakup events (cross section X efficiency) 
plotted as a function of relative energy. The solid and dashed his- 
tograms denote simulated El+Ml and El yields, respectively, (b) 
p-Be coincidence efficiency calculated by Monte Carlo simulations. 
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FIG. 3. The astrophysical Si 7 factors deduced from the present 
experiment plotted as a function of the p- 7 Be relative energy 
(closed circles), in comparison with results from direct and other 
Coulomb-dissociation experiments (closed triangles from |3 , open 
boxes from ^| , open triangles from |J , open crosses from [ 3| , and 
open circles from jjj). The solid curve shows the prediction of 
Bertulani fitted to our data, while the dashed curve shows a fit 
of the theoretical curve of Descouvemont et al. |2o| to the combined 
datasets of Refs. The insert shows the low-energy part of the 

figure, where we have fitted the energy dependence of Jennings et 
al |24j to our data. 
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